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ABSTRACT 

A new class of exact solutions of Einstein’s field equations representing 
anisotropic distribution of matter on pseudo-spheroidal spacetime is obtained. 
The parameters appearing in the model are restricted through physical require¬ 
ments of the model. It is found that the models given in the present work is 
compatible with observational data of a wide variety of compact objects like 4U 
1820-30, PSR J1903+327, 4U 1608-52, Vela X-1, PSR J1614-2230, SMC X-4, 
Cen X-3. A particular model of pulsar PSR J1614-2230 is studied in detail and 
found that it satishes all physical requirements needed for physically acceptable 
model. 

Subject headings: General relativity; Exact solutions; Anisotropy; Relativistic compact 


stars 
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Introduction 


The study of interior solutions of Einstein’s field equations play a significant role in 
predicting the nature of the star in the last stages of evolution. There is an emerging interest 
among researchers to develop mathematical models of superdense stars that are compatible 


with observational data. A number of a rticles appear e d recently in literature tha 


are m 


good agreement with c 

)bservational data 

Pandva et 

ai 

(2014 

); 

Ganeopadhvav et ai 

(2013) 

Hansrai and Maharai 

12006 

); 

Tikekar and Jotania 

(2009 


Baneriee et ai 

(2013 

). The 


spacetime metric representing these models may not have a known 3-space geometry. Many 


models. The spheroidal spacet 

the paraboloidal spacetime us 

ime used bv 

Vaidva and Tike 

kar 1 

1982) 

Tikeka 

I (1 

990) and 

sd by 

Rnch and Skea 119891). 

Tikekar and Jotania 

120071. 

1 ■ 

Sharma and Ratanpal 120131) 1 

^andva et ai 

(2014 

are examples of spacetimes with dehnite 

3-space geometry. The spacetime metric used by 

Tikekar and Thomas 

(1998), 

Thomas et Oi 

12005). Thomas and Ratanpal 

120071 

. Tikekar and .Totania 

2005 

Paul and Chattopadhvav 


fl2010l) spheroidal geometry for its physical three space. It is found that all the above 
spacetimes are useful to describe compact objects like neutron stars and quark stars. 


The matter content of astrophysical object s in ultra 


the form of perfect fluid. Theoretical study of 


li gh de nsitie s may n o 


RudermanI fjl972l ) and 


Canutol 1 19741 ) 


be in 


about more realistic stars suggest that in a very high density regime, matter may not 


be isotropic. Anisotr opy arises due to the exis 


tence of solid stel 


of tvpe-3A sup erfluid 


(ISawver 


Kippenhahn and Weigert 


ar co re or by presence 


(1l99nl ). ISokolovI (1l98nl ). phase transitions 


19721) or pion condensation in a star. Study o 


matter got wide attention after the pioneering work of 


anisotropic distr i bution s of 


Bowers and Liand 1119741) . A 


number of articles appeared in literature related to compact anisotropic stars. The 
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anisotropic model c 
its matter content 


eveloped by 


Maharai 


Gokhroo and Mehra f 


non-uniform density. 


and MarteensI fjl989l) have uniform density for 


Tikekar and Thomas 


99411 gave a m ore realistic model in corp or at ing 


(11998 . 


2005 1. Thomas et ai ( 2005 1 developed 


superdense anisotropic distributions on pseudo-spheroidal spacetime. The spacetime 


used to study the non-adiab atic gravitationa 


pseudo spheroidal geometry. 


Dev and Gleiser 


collapse of anisot ropic distribution has a 


( 12002 . 


2003 


200411 have study the impact of 


anisotropy on the sta bility of stars. Anisotr o pic di st ributions with linear equation of s tate 


have been studied by 


Sharma and Maharai 


Komathirai and IV 


state. 


Su 


models. 


nzu et al 


lara 


(1200711 . 


Thirnkkanesh and Maharai 


(1200811 . 


(120071) studied charged distributions using linear equation of 


(I 2 OI 4 JI us e d line ar eq uation of state for de s cribin g charged quark star 


Feroze and Siddiquil (120111 1 and 


Maharai and Takisal (120121 1 developed models 


of anisotropi c distribu 


quark star^ 


constant. 


i ons u sing quadratic equation of state. In the MIT bag model of 


Paul et al (1201111 have shown that introduction of a nisotro py can affect bag 


Thirnkkanesh and Ragell (1201211 . 


Maharai and Takisal (j2013bll used polytropic 
equation of state f or dev elo ping aniso tropic models. The anisotropic charged star models 


given by 


Malaverl (2013a 


equation of state. 


b 


2014 


Pandva et ai 


and 


Thirnkkanesh and Ragell (120141 1 have Van der Waals 


20141 1 have given aniso t ropic 


with observatio nal data by generalizin g 


model given by 


Sharma and Ratannal ( 201 


th at the models given b y 


by 


Pandva et ai 


Gangopadhvav et al (1201311 . 


compact stars compatible 


Finch and Skeal (1198911 model. The anisotropic 


is a particular case of this model. It is found 


(1201411 agree with the recent observational data given 


In this paper we have studied anisotropic stellar models on pseudo-spheroidal spacetime. 
The geometric parameter R plays the role of radius of the star. The geometric parameter 
K appearing in the model is bounded between two limits to comply with the physical 
requirements. It is found that the anisotropic model developed can accommodate wide 
range of pulsars like 4U 1820-30, PSR J1903-I-327, 4U 1608-52, Vela X-1, PSR J1614-2230, 

















































































































































5 


SMC X-4, Cen X-3. 


In section [2], we have given the spacetime metric and obtained solution of Einstein’s 
held equations. The bounds of model parameter K is obtained in section [3] by imposing the 
physical requirements a physically viable model is expected to satisfy. In section HJ we have 
shown that the anisotropic c ompa ct star models are in agreement with the observational 


data of 


Gangopadhvav et ai fl2013[) and discuss the main results obtained at the end. 


2. Space-time Metric 


We shall take the pseudo-spheroidal spacetime metric for describing the anisotropic 
matter distribution in the form 


ds‘^ = 1 -|- sin^ ddcjP'') , 


1 -I- ^ 

where K, R are geometric para meters and K > 1. The 


( 1 ) 


anisotropic matter distribution (IMaharai and Marteens 


ener gy-momentum tensor for the 


19891) is taken as 


Tij = {p + p) UiUj - ppij + 71, 


IJ, 


( 2 ) 


where p, the proper density, p denotes huid pressure and Ui denotes the unit four-velocity 
of the huid. The anisotropic stress tensor Tiij is given by 


77 , 




= \/ 3 ^ 


CiCj 2 i^UiUj Qij) 


( 3 ) 


Here S = S{r) is the magnitude of the anisotropic stress and c* = (O, —e denotes 

a radial vector. 


We have now the following expressions for radial and transverse pressures 

2S 


Pr = -Tl = p 


V^’ 


( 4 ) 
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and 


P± = -T2 =P- 


S 

7!’ 


( 5 ) 


The difference in radial and transverse pressnres is taken as the measnre of anisotropy and 
is given by 


S = 


Pr - Pi. 

Vs 


The Einstein’s field eqnations 


2 S'KTij , 


are eqnivalent to following set of three eqnations 

1 — e~^ e~^\' 

Svrp = -5-1-. 


,-A 


SttPj, = 


e-V' 


= e ^ ( y + -^ + 


u” _ y. 


4 ' 2r 

Equation (jS jl - ffTOl) can be further modified to the following form, 


A . 2m 

e-^ = 1-, 

r 


2m\ , 2m 

1- ] u = SuprV H- 

j ry>^ 

SttVsS^ = (Syrp + Sirpr) u' + 2 (Svrp^), 


where, 


Using 


m{r) = An v? p{u)du. 


A ^ + 

1 I zi ’ 


in ([8]), we obtain the expression for p in the form 


Svrp = 


K-l S+K^ 




(1 + K^) 


2 \2 • 


( 6 ) 

( 7 ) 

( 8 ) 

( 9 ) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 


(16) 
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The metric potential v can be obtained from equation flT^ once we know the expression 
for -pr- We shall take the expression for radial pressure as 


Snpr = 


K 


1 - 




„2 \ 2 ■ 


(1 + A'f) 

It can be noticed from equation flTTI) that the central pressure is Pr(0) = which is 
directly related to the geometric parameters K and R. This choice of Pr facilitate the 
integration of equation flT^ and obtain in the form 


(17) 


= CR 


k'^-2K+1 




9 N 

\ 2K 


R^J 


\ ^ 


(18) 


where (7 is a constant of integration. 


Differentiating equation (ITT)) with respect to r, we get 


dpr 2r{K - 1) 

OTT—— = -—- 

dr R^ 


It can be noticed from equation (IT^ that 


l + K{2-j^ 


(l + A#)’ 


(19) 


(r = 0) = 0, 

dr 

^ dpr , 2{K — 1) 

(’■ = «) = < O' 


( 20 ) 


( 21 ) 


Further Stt^ < 0 for all values of r in the range 0 < r < 72. Hence the radial pressure i 


IS 


decreasing radially outward and becomes zero at r = R, which is taken as the radius of the 
anisotropic fluid distribution. 


From equation flTB]) . we get 


dp _ 2rK{K-l) 5 + K^ 

dr^ R* (l + K!^f 


(22) 


indicating that the density is a decreasing function of r. 








The spacetime metric ffTlinow takes the explicit form 


ds^ = CR 


X^-2iT+l 


9 \ 

\ 2K 


9 \ 

2 \ 2 


r" \ 




1 

-*- + i?2 


(23) 


-r^ {dO"^ + sin^ Odcj)^') . 


For a physically acceptable relativistic distribution of matter, the interior spacetime 
metric (12^ should continuously match with Schwarzschild exterior metric 


, 9 , 2M\ , . 

ds^ = ( 1-] dt^ 


2M 


-1 

[dO'^ + sin^ Odcj)^) , 


across the boundary r = R. This gives the constants M and C as 

R{K -1) 


and 


M = , 

2{K+iy 

{l + KY 


/ k^- 2K+1 (3K+1\ f5-K\ 

C = Ry ^ ) (1 + 2^—). 


(24) 


(25) 


(26) 


The expression for anisotropy S can be obtained using (fT6|) . flT7|) . fflSj) and (IT^ in ([13]). 
We have 


SttVsS = 


{K - ly 12 + {-6K^ + 16F: - 2) ^ + {-K^ + 3K^ -7K+1) 


r 

B? 


4B4 (1 + (1 + f) 


(27) 


which takes the value zero at r = 0. The expression for Svrpx = Sttp^ — Stt-v/SS', now takes 
the form 


87rp_L = 


{K-l) 4 + (4iF - 12) ^ + (6^2 - 16K - 2) + {K^ - 3K^ + 3K - 1) 




(28) 


3. Physical Analysis 

The anisotropic matter distribution described on the background of pseudo-spheroidal 
spacetime contains two geometric parameters R and K. Since Pr{r = i?) = 0, the parameter 















9 


R represents the radius of the distribution. The bounds on the other parameter K is to be 
determined using the physical plausibility conditions stipulated below: 

(i) p(r), Prir), p±{r) > 0 for 0 < r < i?; 

(ii) p — Pr — 2p_L > 0 for 0 < r < 

(hi) t ^<0for0 <r</?; 

(iv) 0 < ^ < 1; 0 < ^ < 1, for 0 < r < i?; 

(v) The adiabatic index r(r) > for 0 < r < i?. 

The conditions p(r) > 0, > 0, < 0 and ^ < 0 are evidently satisfied in the light of 

equations ffTbjl . fITT)) . fl22|) and flT^ . respectively. 

The condition p_L > 0 imposes a restriction on the value of K, namely, 

K > 2.4641. (29) 


In order to examine the strong energy condition, we evaluate the expression p — pr — 2p_L 
at r = 0 and ai r = R. It is easy to see that 

{p-Pr- 2p±) (r = 0) = 0, (30) 

and {p — Pr — 2p_l) {r = R) > 0, imposes an upper bound for K, namely. 


K < 4.1231. 


(31) 


It has been suggested by 


Canntol (jl974J ) that he velocity of sound should be monotonically 


decreasing for matter distribution with ultra-high densities. This demands that ^ 0- 

The expressions for ^ and are given by 


dp 


dpr _ I+ 2K- Kj^ 
dp ~ K{b + K^) ’ 


(32) 
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It can be noticed from fl3^ that ^ < 1 thronghont the distribntion. 


dp± 

dp 


l + Kj^ 


2K{K - 1) (5 + K^) [2 + + Y^^ + Y^§ + Y,g, + 2K^^] ’ 


(33) 


where, Xi = 8K^+8K-16, X 2 = -4:K^+28K^-20K-4:, X 3 = SK^-AK^-30K^+36K-5, 
X 4 = lOK^ - 3QK^ + IQK^ + 12K - 2 , X 5 = - AK^ + QK^ - AK^ + A and ^1 = 8 ^ + 4, 

Y 2 = 12^2 + IQK + 2, ^3 = 8 K^ + 2AK^ + 8K, Y^ = 2K^ + IQK^ + 12K^, F 5 = AK^ + 8 K^. 


The condition < 1 at r = 0 and r 


R gives the following bonnds on K, viz., 


K > 1.3333, 


(34) 


and 


1<K< 14.7882. (35) 

The expression for adiabatic index T is given by 


T = 


A — — 4- K— 

^ j:j2 W .rv p2 


ij2 


l + 2K-Kj^ 


K (1 - jp) (5 + Kjf) 


(36) 


The necessary condition for the model to represent a relativistic star is that T > 
thronghont the star. T > | at r = 0 imposes a condition on 77, viz.. 


K > -3. 


(37) 


Considering all the relevant ineqnalities, we find that the admissible bonnd for K is given 
by 

2.4641 < 77 < 4.1231. (38) 
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4. Application to Compact Stars and Discussion 


In order to examine the suitability of our model to fit into the observational 
data, we have consi d ered the masses and radii of some well known pulsars given by 


Gangopadhvav et al (120131) . We have considered PSR J1614-2230 whose estimated mass 
and radius are 1.97Mq and 9.69 km. If we set these values in equation (|25|) we get 


K = 3.997 which is well inside the valid range of K. We have further verihed that our 
model is in good agreement with the estimated mass and radii of a number of compact 


stars like Vela X-1, 4U1608-52, PSRJ1903+327, 4U1820-30 SMC X-4 and Cen X-3. The 
value of A, mass, radius and other relevant quantities like pc, Pr, u = ^ and = 0) are 
shown in Table-1. 


Table 1: Estimated physical values based on the observational data 


STAR 

K 

M 

(Me) 

R 

(Km) 

pc 

(MeV fm-3) 

Pr 

(MeV fm-3) 

«(=f) 

II 

o 

4U 1820-30 

3.100 

1.58 

9.1 

2290.97 

277.12 

0.256 

0.465 

PSR J1903+327 

3.176 

1.667 

9.438 

2206.89 

260.52 

0.261 

0.463 

4U 1608-52 

3.458 

1.74 

9.31 

2561.92 

277.50 

0.276 

0.458 

Vela X-1 

3.407 

1.77 

9.56 

2379.27 

261.63 

0.273 

0.459 

PSR J1614-2230 

3.997 

1.97 

9.69 

2883.52 

269.33 

0.300 

0.450 

SMC X-4 

2.514 

1.29 

8.831 

1753.84 

261.06 

0.215 

0.480 

Cen X-3 

2.838 

1.49 

9.178 

1971.21 

260.42 

0.239 

0.470 


In order to have detailed analysis of various physical conditions throughout the star we 
have considered a particular star PSR J1614-2230 having mass M = 1.97Mq and radius 
R = 9.69 km along with the geometric parameter K = 3.997. The variation of density 
and pressure from centre to the boundary of the star is shown graphically in Figure [1] 
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and Figure [21 respectively. It can be seen that density and pressures are monotonically 
decreasing functions of the radial variable r. In Figure El we have shown the variation of 
anisotropy S throughout the star. The anisotropy increases initially and after reaching a 
maximum at r = 2.54, then it starts decreasing till the boundary of the star. The variation 
of square of sound speed and strong energy condition are displayed in Figure [Hand Figure |5l 
respectively. It can be noticed that the square of sound speed is less than 1 and the strong 
energy condition is satisfied throughout the star. 

In Figure [HI we have shown the equation of state for matter distribution in graphical 
form. For a stable relativistic star, the adiabatic index F should be greater than | 
throughout the configuration. We have plotted the graph of F against r in Figure [71 The 
graph clearly indicates that F > | throughout the star. For a physically acceptable model, 
the gravitational redshift, — 1, should be a decreasing function of r. Further the 

central redshift Zc and boundary redshift z^. should be positive and hnite. From Figure [H 
it can be seen that these conditions are satished throughout the star. 


We have studied spherically symmetric anisotropic distributions of matter on 
pseudo-spheroidal spacetime. The model we have developed is in good agree ment with the 


observational data of pulsars recently studied by 


Gangopadhvav et al (l2013l l. The model 


parameters are carefully selected so that the models satisfy all the physical requirements 
throughout the distribution. We have studied a particular model of PSR J1614-2230 and 
have shown that various physical requirements stipulated earlier are satished throughout 
the star. 
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r{Km) 


Fig. 1.— Variation of density against radial variable r. 
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Fig. 2.— Variation of pressnres against radial variable r. 
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r{Km) 


3.— Variation of anisotropy against radial variable r. 
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Fig. 4.— Variation of against radial variable r. 
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r [Km) 


Fig. 5.— Strong energy condition Vs radial variable r. 



Fig. 6.— Variation of pr and p±_ against p. 
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r (Km) 


Fig. 7.— Variation of F against radial variable r. 
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Fig. 8.— Variation of gravitational redshift against radial variable r. 
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